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Abstract 

We consider the pretilt alignment of nematic liquid crystal (LC) on inhomogeneous surface 
patterns comprising patches of homeotropic or homogeneous alignment domains, with azimuthal 
anisotropy assumed in the surface plane. We found that the resultant LC pretilt generally increases 
continuously from the homogeneous limit to the homeotropic limit as the area fraction of the 
homeotropic region increases from to 1. The variations are qualitatively different depending on 
how the distance between patches compares to the extrapolation length of the stronger anchoring 
domain. Our results agree with those previously found in stripped patterns. The present findings 
may provide useful guidelines for designing inhomogeneous alignment surfaces for variable LC 
pretilt control - a subject of much technological interest in recent years. 
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Study of inhomogeneous surfaces for liquid crystal (LC) alignment has witnessed a rapid 
growth in recently years . Experiments showed that these surfaces 

could have unique applications as in the making of multi-stable display devices |7] and the 
control of LC pretilt alignment 0, ^| . This type of alignment surfaces generally comprise 
two kinds of domains favoring different LC orientations. Majority of theoretical treatments 
»d dea, w ith drains _ in a lt e r na ting dipped Q fl fl fl fl « d—ard 
3 S patterns. No systematic calculation has been carried for a patchypattern, which is 
more likely achieved by non-lithographic techniques that are less costly |5|. In this letter, 
we perform an extensive simulation study on the alignment of nematic LC on surfaces con- 
taining patches of homogeneous or homeotropic alignment domains, with one embedded in 
the matrix of the other. 

Figure ^ shows the inhomogeneous alignment surfaces we studied for different area frac- 
tions of the homeotropic region, p. We consider a cubic simulation cell with sides of unity 
filled with nematic LC. The lower (z = 0) surface is the patterned surface as shown in Fig. 
and the upper surface is free. The cell is assumed to repeat in the x and y directions. If 
each type of domain has surface energy Wi(9, 0) (i = 1, 2), the unit-area surface energy, F s , 
would be: 

2 

F a (e(r), 0(r)) = 5{z) £ f^x, y)W l (d(r), <P(r))- (1) 

i=l 

Here, r = (x, y, z) is a position vector in the simulation cell; the LC director field is denoted 
by n(r) = (cos 9 cos 0, cos 9 sin 0, sin 9), where 9(r) and 4>(r) are, respectively, the zenithal 
(measured from the surface instead of the surface normal) and azimuthal angles of the 
LC director; fi(x,y) is one if (x,y) belongs to domain i and zero otherwise so that the 
area fraction of domain i is obtainable by integrating fi(x,y) over the z = unit surface. 
Since our focus is in the LC zenithal (pretilt) alignment produced by the inhomogeneous 
surface, we impose anisotropy in the x-y plane, i.e., both domains prefer the same azimuthal 
orientation along, say, <pQ. We further assume the local preferred zenithal alignment of 
domain 1, # 01 = 88° and that of domian 2, 9 02 = 6° to lift the foreseen degeneracy between the 
4>{r) = 0o and 4>{v) = 0o+18O° states. With these assumptions, the azimuthal orientation of 
the LC is constant in space and equal to 0o- Hence the surface energy depends only on 9(r). 



We use the Rapini-Papoular form of the LC surface energy 



PW,0o) = \W ei sm 2 (9-9 0t ) 



where Wei is the polar anchoring strength of domain i. For the LC elastic energy density, 
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F e , we assume the the Frank-Oseen form: 




(2) 



Here Kn, K22 and K33 are the splay, twist and bend elastic constants, respectively, for 
which the values of the nematic 5CB are used in simulations [8]: K\\ = 7.0 x 10~ 12 N, 
K 22 = 3.5 x 10~ 12 N, K 33 = 1.0 x 10~ n N. If the size of the cubic cell were scaled to A 3 , the 
total energy per unit area is given by 



where F e = F e /Ku, £ ei = K\\/W$i is the extrapolation length of the domain i. Equation (jHJ) 
shows that for the same LC material, the energy density of the system generally increases 
with (1/A), but the critical parameter determining the LC configuration is the ratio X/£ e 
(where £ e is the smaller of £ ei ): If X/£ e <C 1, the elastic energy would dominate the surface 
energy, and vice versa in the opposite limit. For typical alignment surfaces, the polar 
anchoring energy is 10~ 4 — 10 -3 Jm -2 and hence £ e ~ 0.01 — 0.1 /im. We determine the 
equilibrium LC director field by minimizing Eq. (JHJ) with respect to n(r) numerically as in 



We first discuss results obtained in the X/£ e <C 1 limit. Shown in Fig. |2]are the averaged 
pretilt angle of the calculated LC alignment at z = 0, # av (0), plotted versus p for A/£ e2 = 0.01 
and various ratios of Wgi/Wo2 from 0.5 to 5. As seen, all the curves in Fig. |2]bear similar 
shape with # av (0) demonstrating a gentle rise from 6° as p increases from initially, followed 
by an abrupt transition from small to high pretilt whereupon the rise with p becomes gradual 
again. For Wgi/Wg 2 = 1, the transition occurs at p — 1, but shifts to smaller (larger) p 
when Wgi/Wg 2 increases (decreases). 

We plot <jQ y (z), the standard deviation of 8(r) over the x-y plane at z, against z for various 
A/42 horn 0.01 to 1000 Fig. El In these plots, p = 0.125 and Wei = W 92 though we found 
that the results are qualitatively independent of the set values of these parameters. As seen, 
when X/£ e 2 C 1 (= 0.01, 0.1), variations of the LC pretilt in the x-y plane is negligibly small 
for all z, suggesting the director field to be uniform in the present limit. Hence the behaviors 
shown in Fig.|21for LC alignment at z = are essentially the same as those in the bulk. But 




(3) 
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as Xj£ e2 increases towards 1, <Jq V {0) increases notably, and steadily approaches ~ 27° as X/£ e2 
continues to increase towards 1000. We postpone discussion of the noted asymptotic value of 
<Jg y (0) until later. Despite of the large variations in oy^O), all curves display similar decays 
towards within z = A. Least-square fittings reveal that the data of Fig. El can be described 
very well by exponential decays with decay lengths ~ X/2n, consistent with previous results 
found in stripped surface patterns possessing zenithal inhomogeneity 

We return to the data of Fig. EJ In adopting an uniform LC director field, the system 
reduces the elastic energy in the expense of a larger surface energy. The resulting pretilt 
angle should then equal # mm , the pretilt angle that minimizes the surface energy term of 
Eq. (j3J). It is straightforward to show that 9 min satisfies: 

1 = 1 _ Wei sin2(g min - fl 01 ) 
p W 82 sin 2(9 

min — #02) 

In Fig. |21 we plot 9 m i n deduced from Eq. |U versus p for the same set of Wgi/Wg 2 indicated in 
the figure (dashed lines). As seen, the data of 9 min demonstrates excellent agreement with 
those of 6 l av (0) obtained from simulations, verifying our interpretation that the uniform LC 
alignment is due to # min . 

Next we consider the opposite limit where X/£ e ^> 1. In this case, the director field on 
the surface is no longer uniform. But as Fig. El shows, it comprises large surface-induced 
deformations that decays to zero within a distance of ~ X/2n from the surface. Hence the 
homogenized LC alignment in the bulk could be obtained from that at z — A. Shown in 
Fig.HJare plots of 9 ax (X) versus p for X/£ e2 = 1000 and the same set of Wqi/Wq 2 investigated 
in Fig. |21 In the same figure are drawn the corresponding data at z = for comparison. As 
one could see, the average LC pretilts at both z = and z = X are independent of Wgi/Wg 2 . 
It is interesting to note that the z = data demonstrate a simple linear rise with p from 6° 
to 88°. This is because when A/4 3> 1, the surface energy dominates the elastic energy in 
the LC system. One may thus assume that the LC director to exactly follow that of the local 
preferred direction of the anchoring surface, and hence # av (0) should just be the algebraic 
average of local preferred orientation over the surface, i.e., p9 m + (1 — p)9 02 . If p = 0.125 as 
for the cases exemplified in Fig. El 6* av (0) = 16.25°. Then the standard deviation of 9(x,y) 
at z = in the present limit would be 27.1°, in excellent accord with the value shown in 
Fig. El for X/£ 2 = 1000. The slightly larger value of 6* av in the bulk than that at z = is 
attributable to the assumed ^33 being larger than Kn, which favors less bending in expense 
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of more splay deformations of the LC director field in homogenizing into the bulk |2j. 

In this letter, we have considered in detail the pretilt alignment of nematic LC on a 
periodic patchy surface pattern comprising domains favoring either 88° or 6° zenithal align- 
ments. If the pattern period is much smaller than the extrapolation length of the stronger 
anchoring domain, the LC alignment is uniform with the LC pretilt being one that minimizes 
the surface anchoring energy over one pattern period. But if the pattern period is much big- 
ger than the extrapolation length of the stronger domain, the LC alignment at the interface 
copies the inhomogeneity of the surface pattern, which, however, dies out within one pattern 
period from the surface. The resulting homogenized LC zenithal alignment in the bulk is 
approximately given by the average of the local preferred directions of the domains weighted 
by their area fractions in the pattern. Previously Ong et al. (s| created random inhomo- 
geneous surfaces by depositing a discontinuous metal film on a silane (homeotropic agent) 
coated corrugated Si0 2 (homogeneous) substrate, then subjecting the sample to glow dis- 
charge to remove the silane unprotected by the metal islands, before dissolving the metal in 
the end. The average LC pretilt on these surfaces was found to increase monotonically from 
0° to 90° as the average metal film thickness was increased, demonstrating good qualitative 
agreement with our results. Our results also display good consistency with those of Barbero 
et al. |3[ who obtained analytic solutions for LC alignment on stripped patterns favoring al- 
ternatively two dissimilar zenithal orientations (though the two preferred orientations must 
be sufficiently similar to validate their calculations). As most discussions presented here are 
general to LC alignment on inhomogeneous surfaces, it is probably that our results could be 
applicable to more general inhomogeneous patterns with irregular shapes. It is hoped that 
our study would motivate more experiments along the same line of thought. 

We acknowledge support of the Research Grant Council of Hong Kong through project 
no. HKUST6115/03E. 
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FIG. 1: Patterns used to simulate the inhomogeneous surface. Filled squares stand for the 
homeotropic domains; unfilled squares stand for the homogeneous domains. Periodic boundary 
conditions are imposed in the x and y directions. 

FIG. 2: Simulation results of average pretilt angle 8 !xv {z = 0) versus p for X/£ e 2 = 0.01 and various 
Wq\/Wq2 = 0.5 (crosses), 1.0 (triangles), 2.0 (diamonds) and 5.0 (squares) are shown. Dashed lines 
are the simulated 9 m \ n versus p according to Eq. (j3J) for the same set of Wq\/Wq2- 

FIG. 3: Profiles of the standard deviation of the LC pretilt along z for different ratios of \/(. e 2 
from 0.01 to 1000. In these data, p = 0.125 and Wei = W e2 . 

FIG. 4: The average LC pretilt at z = X, which is essentially the bulk value, versus p for {\jt e 2 = 
1000) for four different W$i/W$2 as indicated. The data for average pretilt at z = are also drawn 
(solid line). 
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Fig. [U by Wan et al. 
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Fig. El by Wan et al. 
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Fig. El by Wan et al. 
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Fig. H by Wan et al. 
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